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Abstract 



We solve 0{n, 1) nonlinear vector model on Bethe lattice and show that 
it exhibits a transition from ordered to disordered state for < n < 1. If 
the replica limit n — > is taken carefully, the model is shown to reduce 
to the corresponding super symmetric model. The latter was introduced by 
Zirnbauer as a toy model for the Anderson localization transition. We argue 
thus that the non-compact replica models describe correctly the Anderson 
transition features. This should be contrasted to their failure in the case of 
the level correlation problem. 
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I. INTRODUCTION 



In the study of systems with quenched randomness historically the first way of treating 
the disorder averages was the so called replica method of Edwards and Anderson In this 
method one introduces n "replicas" of the original system and calculates annealed averages 

for this replicated system. Then the use of identity 

- , Z^- 1 

log Z = lim 

allows to recover the properties of the original system. 

For the problem of Anderson localization replica fields describing electrons may be chosen 
to be either fermionic (Grassmann) or bosonic. In either case one can construct effective 
field theoretic description in terms of non-linear cr-model for matrix field which belongs to 
some compact manifold for fermionic replicas and some non-compact manifold for bosonic 
ones (see, e.g. p). 

Another way to treat the disorder is the supersymmetry method of Efetov [Q. In this 
method one introduces both fermionic and bosonic degrees of freedom and the resulting 
(j-model field is of the supermatrix structure. 

It was realized some time ago that the replica method suffers from serious drawbacks. 
Verbaarschot and Zirnbauer showed explicitly [Q that the replica method fails to give correct 
non-perturbative result for a problem of energy level correlation which is equivalent to 
zero-dimensional cr-model, while the supersymmetry method works nicely. Since then the 
prevailing opinion in the literature on the Anderson localization problem was that the replica 
method may be at best considered as perturbative tool not being able to describe properties 
of disordered (localized) phase. 

The goal of our paper is to reconsider the correspondence between supersymmetric mod- 
els and non-compact models in the replica limit. For this purpose we analyze the solution 
of two "toy" models on Bethe lattice. The first one is one of the simplest models with non- 
compact symmetry, namely 0(n, l)/0(?2)-vector model. The other is its supersymmetric 
counterpart, the so called "hyperbolic superplane" (see 
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The paper is organized as follows. In Sec. H we set up the general description of 
properties of the 0{n, l)-model in terms of a single "distribution function of local order 
parametr" P{6). Doing that in the spirit of dimensional continuation we consider parameter 
n to be any real number large enough to ensure convergence of integrals. In Sec. |T| we 
similarly consider supersymmetric version of our model. Sec. |rV| is the central section of 
the paper. There we discuss the replica limit n — and show that if we take it carefully all 
the results for the 0{n, 1) model exactly reproduce results of the supersymmetric treatment. 
In the following sections we proceed to solve the 0{n, 1) model. Our analysis follows very 
closely that of previous papers devoted to the problem of Anderson localization on Bethe 
lattice and therefore we omit many details. We find that 0{n,l) model exhibits two 
phases with a phase transition between them for any < n < 1. We obtain critical behaviour 
of different correlators near this transition and show that it is exaclty the same as exhibited 



by the supersymmetric model. Finally, Sec. VIII contains discussion of our results. 



II. 0(n, l)-MODEL: GENERAL EQUATIONS 



We start with Hamiltonian 



H = J ^ iii ■ Uj + H^ai. 

<ij> i 

Here i and j refer to the sites of Bethe lattice with coordination number m+1 and n = (a, tt) 
is a (n + l)-component vector sweeping the hyperboloid if"'^, defined by the equation 
= cr^ — 77^ = 1. This hyperboloid is the symmetric space, associated with 0{n, 1) group: 



jjn,i _ Q(ji^ 1)1 0{n). We parametrize n as follows a = vT+Tr^ = cosh 6', < 6 < oo. 



TT 



sinh 6 cos 0i 
sinh 6 sin 0i cos 02 

sinh 6 sin 0i sin 02 • • • sin ( 



)i, . . . ,0„_2 e [0,7r], 0„_ie[O,27r) 
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With this parametrization the scalar product is rij ■ nj = aiCj — tTj ■ ttj > cosh(^.j — 6j) > 1 
and, therefore, the hamihonian Ti is bounded from below only for J, H > 0. The 0{n, 1)- 
invariant measure on if "'^ is 



where a is a normalization constant to be fixed later. 

Now we introduce "distribution function of the local order parameter" P(n) in the usual 
manner. Namely, we cut one of the m + 1 branches coming from site n and integrate the part 
of the Boltzmann weight exp(— 7i) over this branch. The resulting function P(n) satisfies 
the integral equation 



where we introduced the following notation: L(n, n') = e~'^" "', -D(n') = e'^'^' . Knowl- 
edge of the function P(n) allows us to calculate the partition function Z = JdnD{n)P"^^^{n), 
one-site averages {A{n)) = Z^^ JdnA{n)D{n)P"^^^{n), and "weighted" correlators 



where the factor N{r) = (m + l)m'''~^ counts the number of sites located at the distance r 
from a given site (without this factor all the correlators exponentially decay because of the 
geometry of the Bethe lattice). 

The constant a in the definition of the measure dn can be chosen arbitrarily. It is easy 
to see that rescaling of the measure changes overall normalization of -P(n) and Z but does 
not affect either one-site averages or correlators. This allows us to choose a convenient 
normalization for P(n) as follows. Note that when H = the Eq. (|l]) admits constant 
solution. Then we require that this solution be simply P(n) = 1 or, equivalently, that 



dn = a d9 sinh" ^9 d(j)i sin" d(j)2 sin" ^4>2 - ■ ■ d(j)n~i 




(1) 





M(n,n') = L(n, n')D{n')P"'-\n'), 




(3) 
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This fixes a = [0'" {2K^^_,y2{J))-\ 

Magnetic field H breaks the 0{n, 1) symmetry down to 0{n). Then the function P may 
depend only on cr or, equivalently, on 6. This allows us to perform angular integrations in 
Eq. dg), yielding 

roo 

p{9)= de'LL{e,e')D{e')P'^{e'), u) 

Jo 

r n/\ ^ I ^^i^'^ \ ^- J cosh 6» cosh 9' 

X {271 J sinh 6 sinh 9')^Ir_i{J sinh 6 sinh 6'), 
D{d') = e-^-^t^^'. (5) 

Similar integration may be done in expressions for partition function Z = 
aSn~iI^desmh''-^eD{e)P"'+^{e), where Sn-i = 27r"/Vr(n/2) is the volume of the sphere 
5"^^, and correlators. In particular, upon averaging of n only a-component survives giving 
the "order parameter" 

(a) = (cosh^) = — — / desmh''-^ecosheD(e)P"'+\e). 

Z Jo 

For invariant correlator (ng ■ n,.) = (aQCXr) — (ttq ■ tt,.) the angular integrals give different 
kernels for longitudinal G'^i^r) = (cosh 6*0 cosh 6'r)«, and transverse Gfj{r) = {7iQiHrj)w parts: 

G^(r) = /"rf^osinh"-^0ocosh0o^(^o)^"(^o) 
Z Jo 

X (jll^de,ML{e,_^, e,)^ P{9r) cosh^, (6) 

Tl/j Jo 

X (f[l°^deiMT{ei.i,ei)^p{er)smher, (7) 

where 

MLi9, 9') = LLi9, 9')D{9')P"'-\9'), Mt{9, 9') = Lt{9, 9')D{9')P"'-\9'), 
T (a Q/\ 1 /sinh6''\'" ^- j cosh e cosh e' 

X (27r J sinh 9 sinh 9')^!^ (J sinh 9 sinh 9'). 



Equations (||), (|^) and may be written symbolically as {A{no)B{nr))w = 
Z^^N{r){A{nQ)\M^\B{nr)), where M represents an integral operator with one of the kernels 
M, Ml, or Mt- Introducing complete set | (f)\) of eigenfunctions of M: M \ (f)\) = A \ 
(f)js) we can rewrite correlators as 

We chose operators M to be non-symmetric, which means that left and right eigenfunction 
are different. 

III. SUPERSYMMETRIC VERSION: HYPERBOLIC SUPERPLANE 

In this section we consider the supersymmetric version of 0{n, 1) model, namely, a non- 
linear model with field taking values on the so called hyperbolic superplane. This object is 
constructed as follows (see We consider a set of 5-component vectors 

^ = (a,7ri,7r2,^, -0 

where the first three components are commuting, whereas the last two are Grassmannians 
(we use the adjoint of the second kind, see llO| for a review of supermathematics). Next we 



consider the group G of linear transformations in the space of vectors ip which preserve the 
"length" ||'0p= 0"^ — ttJ — 7r| — 2^^. Let K be the subgroup of G which separately preserves 
0"^ and TT^ + 7r| + 2,^,^. Then the coset space G/K is isomorphic to the space of vectors ip of 
unit length ||'?/'||= 1. This is the hyperbolic superplane. 
We will use the following parametrization of G/K: 

I ~2 ~ 

7Ti = sinh 6 cos 0, = sinh 6 sin 0, a = yl + sinh 6 + 2,^,^ = cosh 6 + 



cosh 6' 

In this parametrization the G-invariant measure on G/K is 



I _ / 

dip = a—d7rid7r2d^d^ = a 1 ^ dO sinh Odcpd^d^. 

a \ cosh 6. 



The Hamiltonian in this case is taken to be 



7^ = jjZM + ^E^- 

<ij> i 

We again choose the constant a in the definition of dip such that / dip' exp{— Ji/jiIj') = 1. This 
gives a = e"^/27r. Proceding hke in section |I| we introduce function P{ip) (by symmetry it 
actually depends only on a) which satisfies the equation 

P(a) = jdiP'e-^^^'e-"'''P'^{a'). 

Expanding both the left-hand side and the right-hand side in powers of Grassmann variables 
and integrating them out we get from the last equation 

/■oo 

p(^) = e^(^-=°^^^)e-^p'"(o) + / de'LLo{e,e')D{e')P"'{e'), (9) 

Jo 

LLo{e,e') = e^Jsinh^e-^'=°'^^^°'^^7i(Jsinh^sinh^') (10) 

The first term in (P) is the boundary term resulting from integration by parts. 

If we put 6' = in Eq. (|) we obtain P(0) = e~^P'^{0) which means that P(0) = 
or P(0) = 0. To have P{9) = 1 as a solution for if = we have to choose P{0) = e^-i . 

We can also perform Grassmann integrations in fromulae for partition function, one-site 
averages and "longitudinal" correlators: Z = e'^e'^^, {A{a)) = A{0), (A((To)-B(o"r))«) = 
N (r) A{0) B (0) . In particular we have 

(cosh^) = l, G^(r) = (aoa,)^ = iV(r), (r) = 0, (11) 

where subscript c refers to connected correlator. For transverse correlator Gfj{r) = {'Koi'Krj)w 
we obtain after some calculation 

Gl{r) = 6^^^ j^de,D{d,)P"\d,) 

X (i{J^ de,LTo{e,.^,e;)D{9;)P^^-\e,)^ p(0,)sinh^,,, (12) 

LToiO,e') = e*'Jsinh^e--^^°^^^=°^^^7o(Jsinh^sinh^')- 
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IV. REPLICA LIMIT 



Now we study what happens with the equations of Sec. |T| for the 0{n, l)-model in the 
rephca hmit when n ^ 0, namely whether they reduce to those of Sec. [TT^or not. 

First of all, if we simply set n = in Eq. (§ we get P{e) = f dO'LLoie, 9')D{9')P''^+^{e') 
with kernel (PUf), which differs from Eq. by the absence of the boundary term. From that 
we could conclude, in particular, that -P(O) = in the replica limit and, therefore, this limit 
gives the incorrect answers. However, this simple recipe is wrong. To see that, let us set 
6* = in Eq. (|p before taking the replica limit. Using small- 2; expansion Iu{z) ^ r{u+i} (f) 
valid for u 7^ —1, —2, we get 

POO 

P(0)=aSn-i dO sinhJ"-^ 66-^"°'^^ D(e)P"'(e). (13) 
Jo 

If we assume that P(0) 7^ then integral in the last equation is divergent at the lower limit 
if we set n = 0. At the same time a Sn-i = for n = 0, so the expression ([T3| ) is of the type 
• 00 in the replica limit and should be studied in a more careful way. For this purpose, let 
us consider the following identity: 

roc roc roo 

Sn-1 / (i^sinh^-i 9/(9) = /(0)5„_i / ci^sinh"-^ 9 + S^-i / (i^sinh"'^ 9{f{9) - /(O)). 
Jo Jo Jo 

Here the first term is finite for < n < 1 and gives 7r*^"^^^/^r (''^) /(O)) whereas in the 
second term the difference f{9) — /(O) makes the integral convergent even for n = 0. Now 
we can safely take the replica limit in which the second term disappears, and we get 

roo 

lim5„_i / d9smh.''-^9f{9) = f{0). (14) 
Therefore, the replica limit of Eq. (|13D is simply P(0) = e~^P"^(0) which is exaclty the 



result for P(0) from Sec. pi 



Now we can perform similar trick for arbitrary 9: 

roo roo roo 

/ d9'LL{9,9')f{9') = f{Q) d9'LL{9,9')+ d9'LL{9,9'){f{9') - f{0)). 
Jo Jo Jo 
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Before we take the rephca hmit, the integral in the first term here equals 1 due to normal- 
ization of kernel L(n, n'), Eq. (|^). In the second term we can safely take replica limit as 
before. The kernel there becomes Llq{6, 9') of Eq. (p^O]). After that we split the second term 
again into two pieces to get 

hm / de'LL{e,e')f{9') = f{Q)(i- d9'LUd.o')]+ de'LUd.e')f{9'). 

n^OjQ \ Jo J Jo 

The integral in the first term can be done (see, e.g. and we get as the result 

poo fOO 

hm / de'LL{e,e')f{e') = /(o)e-^(i-™^'^^) + / de'LLo{e,e')f{e'). (15) 

n^QJo Jo 

Then in the replica limit Eq. (^) becomes exactly the Eq. of Sec. |T|! 

Alternative way of getting Eq. (^) is to separate the two contributions to the modified 
Bessel function which enters the integral kernel Ll{0,6') using the recursion relation 

n 1 /z\^^^ 



The crucial point is that the first term here cannot be neglected, although it has a vanishing 
coefficient in the limit n — 0, since the corresponding integral over 6' will diverge in this 
limit. Thus, in full analogy with the treatment of Eq. (0), this singular contribution should 
be first evaluated at finite n, and only then the limit n — can be taken. This again gives 
the Eq. (|15|). 

The use of prescriptions (0) and (|1^) allows us to show that equations for all the quan- 
tities of interest from the Sec. |I| in the replica limit reduce to those of the supersymmetry 



method of Sec. |ITT[ This is the main result of our analysis. In the remaining sections we 
find that our model exhibits ordered and disordered phases for < n < 1. We find the 
transition point between them and solve for the critical behaviour of correlators near this 
transition. We explicitly show then that in the replica limit this behaviour is identical to 
the one found in Refs. J^,^ for supersymmetric cx-model and in [^] for the Anderson model 
on Bethe lattice. 
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V. ORDERED PHASE AND TRANSITION POINT 



For J ^ 1 we expect an ordered state with spontaneuosly broken symmetry, where all 
the n's are slightly fluctuating around a-direction. The transverse components tt are small 
and we can expand in them: a ~ 1 + ^tt^ + ■ ■ ■. Then the kernel of the Eq. (|l]) becomes 
Gaussian and, therefore, the equation admits (for any value of n) a solution which is also 
Gaussian: 

P(n) = ce-5°^', (16) 

where a = ^ [j{m - 1) - H + ^( J(m - 1) - H)^ + AmJH) {a = for H = 0) and c 

is some normalization constant. With this form of -P(n) we perform Gaussian integrals to 
flnd for small magnetic fleld 

G^{r) oc ^^(r) = exp (-^;;^^-) , (17) 

where "spin stiffness" ps = J. So longitudinal correlator is massive, while transverse modes 
are Goldstones with long-ranged correlations for H = 0. 

Non-compact nature of the 0{n, 1) model makes it rather unusual from the traditional 
point of view. For example, mean fleld theory for this model does not give any transition at 
all: the system is always ordered. If we expect to have disordered phase for J <^ 1, then in 
this phase the order parameter will be large for small magnetic fleld, diverging when H ^ 0. 

The symmetry breaking factor D{6), Eq. (§), becomes signiflcant for cosh^^ ~ H^^ or 
~ ln-|:, therefore it is convenient to change variables to t = ln(if cosh^) in Eq. (^. 
After such a change the argument of Bessel function becomes large for small H, and, using 
asymptotic formula (2tczY^'^I^(z) we arrive in the limit H ^ at 

/oo 
dt'L{t-t')D{t')P^{t'), (18) 
-oo 

= ^ — exp — — t- Jcosht , D{t) =exp(-e*). 

This equation was derived under the assumption that the solution P is a function of variable 
H cosh 6 for small H. Such a solution becomes P = 1 in the limit H = 0. Solution of the 
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"ordered type" (|T6D is not contained in the Eq. (0) or, rather, it corresponds to trivial 
solution P = 0. 

The qualitative behaviour of the solution of the Eq. (|18[) should be as follows. For large 
positive t [t ^ 0) the function P{t) exponentially goes to zero because of the symmetry 
breaking term D{t'). For large negative t (t ^ 0) the main contribution to the integral in 
Eq. ( |18|) comes from t' ^ because the kernel L{t — t') is sharply peaked at t = t'. But 
in this region D(t') ~ 1 and the Eq. (|l^) admits solution P = 1. In the region t ~ there 
should be a kink connecting the two asymptotic solutions. The precise form of the solution 
can be found numerically by iterations starting with P = 1. But such iterative procedure is 
convergent to the solution of described type only for small enough J. If J > J^, where Jc 
is the critical coupling, this solution becomes unstable, the kink at t ~ starts to move to 
the negative t until the solution P{d) takes the form characteristic of the ordered phase. In 
fact, Eq. (|l^) is very similar to the one studied by Zirnbauer [Q, so we follow very closely 
the stability analysis of this paper. 

For t ^ we drop D(t') and linearize the simplified equation around constant solution: 
P{t) = 1 - 6P{t) with 

/oo 
dt'L{t - t')5P{t') (19) 
-oo 

Because of the translational invariance of this equation it admits exponential solutions: 
JZcdt'L{t - t')e^'' = A{u)e''\ where 

A(-) = (20) 

is an eigenvalue of the integral operator with the kernel L(t — t'). 

Let us make a mathematical remark. In terms of vector n Eq. (|1^) corresponds to 
5P(n) = m Jdn'L[n, n')6P{n'), i.e., 6P{n) is an eigenfunction of an 0{n, l)-invariant integral 
operator L with kernel L(n, n'). The space of such functions is spanned by the so called 
spherical functions of group 0{n,l) (see, e. g. iT^). In our case of 0(n)-symmetric P{0) 
these are the zonal spherical functions 
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M0) = (sinhe)i-"/2p;;„/;+\(cosh^), (21) 

where Py{z) is the Legendre function. The funciton il)v{9) is an eigenfunction of L with an 
eigenvalue A(z/) given by Eq. pO|). Some properties of these functions are presented in the 
Appendix. 

For the problem of stabihty of the asymptotic solution P = 1 the relevant values of u 
are real positive numbers. Indeed, the natural perturbation SP{t) induced by the symmetry 
breaking term D{t) ^ 1 — e* in the region t <^ is SP{t) = e*^* with u = 1. Analysis similar to 
that of 1^ shows that if mAmj„(z/) > 1 then the solution P = 1 is unstable and collapses to the 
trivial solution P = (ordered phase). This happens for any value of coupling constant J for 
n > 1 because in this case K(^n-i)/2+u{J) > K(^n-i)/2{J) and mAmini^) = rnA{0) = m > 1. 
On the other hand, for < n < 1, K(^n-i)/2{J) = -^|(n-i)/2| (<^)5 and + u may be 
smaller than | | . In this case there is a transition at the value of the coupling constant Jc 
determined by the equation 



m- 



1. (22) 



K\(n-l)/2\{Jc) 

Solution of this equation for m = 2 is shown on Fig. 1. 

VI. CORRELATORS IN THE DISORDERED PHASE 

From now on we restrict our attention to the case < n < 1. We will also assume that 
magnetic field is very small: H <^1. In this case we can approximate P(0) by 1. Then for 
J < Jf. function 



p{e) 



1, ^<ln| 
0, ^>ln^ 



with a kink connecting these two regions near 6' = ln-|:. Typical solution of Eq. (^) of this 
type for n = 0, J = 10^^, H = 10"^° is shown on Fig. 2. Approximating this solution by a 
step function we estimate partition function as 
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Jo 



which is not singular as — > 0, so we'll calculate it for = 0: 

Z = a Sn^i Hde sinh"-^ e = a7r("-i)/2p fl^] . (23) 



2 

The finitness of this quantity reflects the fact that the total volume of the hyperboloid H"^'^ 
is finite for < n < 1. We also get the "order parameter" 



{a) ^ / de sinh"-^ e cosh 9 ^ (24) 

Z Jo nZ 
So (cr) diverges when H in the disordered phase, unless n = 0, in which case (cr) 
becomes non-critical (similar to the density of states in Anderson transition). 

To obtain the expression for correlators we again perform the change of variables 
if cosh 6* = e* in Eqs. (H), (0). To the leading order in both and Gjj become 
the same (up to a factor Sij/n): 

aSn^iN{r) 



G"(r) 



dtoe^''>D{to)P"'{to) n / dUML{U-i,U) P{U)e' 

-oo \ J— oo / 



z 

Range of integration for operator with kernel ML{t,t') is infinite, which means that the 
spectrum of eigenvalues is continuous. In this case Eq. (|^) for the correlator G{r) is 

G^(r) = (^ + iy"-i ff-n-i rdXW{X)A\X){mAxy, 
mZ Jo 

where A'^{X) = {e^^\4'xi't)){4'\{'t)\^^) ^-^^d W{X) comes from the normalization of eigenfunc- 
tions. Right and left eigenfunctions (t>^x(t) in the limit H ^ should approach ipi^+i\{^) 
and sinh"~^ 0ijji^_^_^)^{6), or asymptotically for t ^ 

fx\t) oc sin (At + 5(A)) (25) 

A 

with eigenvalue = K^^^-!f/l{j) ■ suggested by Eqs. (^) and (|A5|) the normalization 

of 4>x\t) should be 

dtcp[mi,{t)-^^^~^'^ 



W{X) 
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where W{X) oc for small A. Small A behaviour of "phase shift" (5(A) in Eq. p5|) is found 



matching the asymptotic behaviour p5| ) to (f)\{t) = for t > 0, which is again the effect 
of the term D{t). Such matching gives that 6{X) should be at least linear in A for small A. 
Using this we can show that A{X) ~ 0(1) for small A. Also expanding mAx m Jc ~ J and 
in A we find mkx ~ 1 — ao(t/c ~ J) — a2A^. Combining all the above results we arrive at 

JO 

Gj.(r) oc 5.,H-^~\-'/'e-^'^, 

where ^ ~ (^c — J)^^- 

In the replica limit these equations reduce to G^{r) = and 

which has the same form as the density-density correlator in the localized regime in Refs. 

VII. CORRELATORS IN THE ORDERED PHASE CLOSE TO THE 

TRANSITION 

In the case J — Jc Jc (just above the transition) we expect spontaneous symmetry 
breaking which in terms of function P{6) takes place at some large scale A divergent at 
J = Jc such that 



Pie) 



1, < ln2A 
0, > ln2A. 

Similarly to the previous section we find that partition function is nonsingular as A ^ cxd, 
so we take it to be equal to its value at Eq. (ESI). Again, like in Eq. (121) we find 



(a) ^ (26) 
nZ 

The long distance behaviour of the correlators in this phase is determined by the eigen- 
values of the operator with kernel M(n, n'), see Eq. (|]). Because of the form of the function 
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P{0) the integration range in 6 for this operator is finite, and the spectrum of its eigenvalues 
is discrete. It is easy to see that for H = the largest eigenvalue of this operator is 1/m. 
Indeed, for Q G 0{n, 1) using invariance of the kernel L(n, n') and the measure dn we get 
P{Qn) = Jdn'L{n, n')P™(Qn'). Taking Q to be infinitesimally close to unity and expanding 
we obtain 

^^TT = m [dn'L{n,n')P"'-\n')^^^7z'. (27) 
da J da' 

The functions 1/,) = tt^ are the Goldstone modes associated with the symmetry break- 
ing. Integrating out angular variables in (|27|) we have 

In the asymptotic region K 6*, 6^' < In 2A the kernel Lt(6', 9') ^ L{9 - 9') of Eq. (|T8]), Eq. 
(PSP becomes the same as Eq. ([T^) which means that P'{9) has the asymptotic behaviour 



1 — n I 



where A satisfies -^'a(J) _ j_ Expanding this in small A and small J — X we find 

A~(J-J,)i/l (30) 



The fast decrease of ^{9^ near = In 2/1 chooses the value of A such that function (pQ]) 
has the first node a this point: A = j^^. Combining this with (pO]) we find 

A ~ exp (const (J - J^)"^^^) • (31) 

The value of the constant C in Eq. (^) may be found as follows. dP/d9 is re- 
lated to the function 6P{9) of Eq. ( [19|) in an obvious manner: ^ = Using 
this and writing 6P{9) ^ Ci^^^^ sin(A6' + 5(A)) we get 5(A) = tan^^ (^^j ^ ^ and 
C = —Cwjx^ + (^^)~ ~ ^C*!!^^!. But function (5P should approach the value 1 near 
9 = \n2A. That gives Ci ^ |^|(2A)("~i)/2 and 

C^-(^!i_l)'(2A)("-i)/2. (32) 
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Now we evaluate the largest eigenvalue Kmax of operator with kernel M (n, n') for small 
H using the first order perturbation theory: K^ax = {fi\M\fi)/{fi\fi) (no summation!). 
Since kernel M(n, n') is non-symmetric, differs from by the factor P™^^(n). Then 
we have 

2 ^ / , \ 2 



Using Eqs. (|^) and (|32D we estimate to be 

^ ' ' n \^ Jo n 
In the presence of H Eq. (^) is replaced by 

dP(n)^ = Jdn'L{n,n')^ (p™(n')e-^'^') tt'. 



Using this we perform integrations by parts and keep only linear terms in H to find 



m rn?{m + 1) 

Then for the maximal eigenvalue we have mKmax = 1 — H/{{m — l)ps), where the "spin 
stiffness" 

Ps = -T — 7^r\- ^ — A W - Jc) ' exp ( -const J - Jc M • 33 

For the final evaluation of the correlator (A(no)-B(n.r))^ we have to calculate A^(A) = 
I]iLi(^|/j)(/j|-B)/(/i|/i). For longitudinal correlator A(n) = -B(n) = cosh 6' in which case 

cosh^) = /<inP'"(n)a^^7ri = and ^^(A) = 0. This means that the decay of G^(r) is 
governed by smaller eigenvalues and, therefore, G^ir) is massive even for H = 0, which is 
to be compared with Eq. ([T7|). On the other hand, for Goldstone modes we have (/i|vrj) = 
/(jnP"'(n) '^^J'^^ TTiTTj = —Sij^:^ and for the transverse correlator we finally have 



Note that again, as before, in the replica limit our equations (0), and ( PH ) exactly 
correspond to the results of Refs. • 
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VIII. DISCUSSION 



In this paper, we have solved the non-compact 0{n, 1) model on the Bethe lattice for 
arbitrary n. The analytical continuation procedure allows us to consider n as being an arbi- 
trary positive number. We find that for n > 1 the symmetry of the model is always broken, 
so that the system is ordered and the order parameter has a finite value. This is in agree- 
ment with Ref . [0 , where the system was shown to be ordered for n = 1 and n —>■ oo above 
two dimensions. (Note that the Bethe lattice is effectively infinite dimensional). However, 
for < n < 1 we find a transition from the ordered to disordered phase, when the coupling 
strength J decreases. In the replica limit n — 0, our solution reproduces exactly that for the 
supersymmetric version of the model. The latter is nothing but the hyperbolic superplane 
introduced by Zirnbauer as a toy model for the Anderson localization transition. In the 
whole region < n < 1 the qualitative picture of the transition and the critical behavior 
are analogous to those of the supersymmetric model of Zirnbauer, which shows in turn all 
the essential features of the Anderson transition on Bethe lattice studied in Refs. • The 
success of the replica trick may seem surprising, since we know [^] that it gives wrong result 
in the case of the level correlation problem. The crucial difference is that near the Anderson 
transition only the non-compact sector of the supersymmetric cr-models is essential, whereas 
for the level correlation problem both compact and non-compact sectors are equally impor- 
tant. The similar reason explains the agreement between the replica trick renormalization 
group calculation of asymptotics of various distributions and the recent study of these 
asymptotics via the supersymmetry method ||T5[| . 

When the manuscript was in preparation, we learnt about the work by T. Dupre [T^ who 
studied the supersymmetric model of Zirnbauer in 3D and found the critical behavior 
analogous to that expected for the Anderson localization transition. 
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APPENDIX: 



The functions (^) have the following asymptotic behavior 

M0 » 1) 



2n/2-l / Y(u 



ihizi] r(—i/ — 'i^) 

2 ' ^ue ^ ^ \ ^ 2 ' -(v+n-l) 



\T{v + n-l) 



(Al) 



Among the functions (|21|) there is a special subset with v = + , A > 0. These 



functions are real and form a continuous basis in the space (^[1, oo), ti^sinh"' ^ Q^: 

roo 1 

(i^sinh--^ ^^i^+a(^)^V+^A'(^) = - A'), 

POO 

/ d\n{\)iji^^.^{e)i;i^^^^{e') = sinh^-" 05(9 - 9'), 

Jo 2 2 



(A2) 



(A3) 



/x(A) 



r(^ + a) 



(A4) 



For n ^ 1 and for small values of A, the asymptotics ( |A1| ) and the spectral measure ( |A4| ) 
take the form 

2"/2 „o/n-l 



A 



sinA^, /x(A)^r' 



A^ 



(A5) 
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